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NOMENCLATURE 


C,  constant.  Equation  13 

D,  constant,  Equation  15 

E,  modulus  of  elasticity,  Equation  5 
f E],  matrix.  Equation  27 

[F],  matrix.  Equation  27 
fC],  matrix.  Equation  27 
Kj  f constant.  Equation  (A-l) 

K 2,  constant,  Equation  (A-l) 

M,  constant,  Equation  14 

N,  constant,  Equation  14 

P,  applied  longitudinal  load,  Equation  11 
T,  temperature  rise,  Equation  5 

a,  inner  case  radius.  Equation  11 

b,  outer  case  radius.  Equation  11 

c,  inner  chamber  radius.  Equation  24 

d,  outer  chamber  radius,  Equation  24 
£ti  total  radial  strain,  Equation  1 

fc'i  total  axial  strain.  Equation  4 

t. ft  total  tangential  strain,  Equation  2 

• A 

equivalent  plastic  strain,  Equation  8 
yield  stress  in  shear.  Equation  7 
p,  propellant  pressure.  Equation  20 

r,  radial  coordinate,  Equation  1 

u,  radial  displacement,  Equation  3 

s,  axial  coordinate,  Equation  4 


2 


NOMENCLATURE  (Cont'd) 


*X  coefficient  of  thermal  expansion.  Equation  5 

^)t  constant.  Equation  16 

constant.  Equation  16 
axial  stress.  Equation  5 
OV  radial  stress,  Equation  1 
^0  tangential  stress.  Equation  1 
0 tangential  coordinate.  Equation  1 

It  constant,  Equation  8 

P Poisson's  ratio,  Equation  5 

Superscripts 


elastic.  Equation  4 
plastic.  Equation  4 
chamber.  Equation  A-l 


INTRODUCTION 


The  objectives  of  the  Small  Arms  Project  1J562604A607  are  to 
provide  exploratory  development  of  new  or  improved  munitions  com- 
ponents, and  ammunition  simulation.  Task  25  of  this  project  deals 
with  the  mathematical  modeling  of  ammunition,  and  part  of  the  work 
performed  under  that  task  is  reported  here. 

This  study  was  an  extension  of  the  modeling  of  case  extraction 
described  in  an  earlier  report*  by  the  author.  In  that  report  an 
analytical  model  for  determining  the  force  required  to  extract  a spent 
cartridge  case  from  a conventional  weapon  was  presented.  It  was 
recommended  in  that  report  that  additional  refinement  of  the  mechanics 
in  the  model  be  performed.  Part  of  this  additional  refinement  is  pre- 
sented herein. 

In  Reference  1 the  cartridge  case  was  modeled  as  a thin  membrane 
shell  being  deformed  and  stressed,  initially,  by  the  mechanical  pres- 
sure from  burning  propellant  gases,  and  subsequently  by  interaction 
with  the  inner  chamber  wall.  The  chief  purpose  of  the  model  was  to 
calculate  the  state  of  residual  radial  stress  in  the  case  at  extraction. 
This  calculated  residual  radial  stress  at  the  outside  diamete  * of  the  case 
wall,  when  multiplied  by  the  area  on  which  it  acts  times  t coefficient 
of  friction,  and  less  the  blowback  thrust,  was  essentially  the  force  of 
extraction.  (See  Equation  6 of  Reference  1).  Those  factors  which 
were  neglected  or  over-simplified  in  that  model  and  are  corrected 
here  are: 

1.  The  ability  to  include  thermal  (steady  state  and  transient) 
strains. 


2.  The  compatibility  of  the  radial  and  hoop  strains. 

3.  The  use  of  simplified  equilibrium  and  constitutive  equations  in 
which  the  state  of  street  was  independent  of  the  elastic  and  plastic 
material  properties  of  the  case. 


1 

P.  Gordon,  "Analytic  Study  of  Extraction  in  the  M16  Weapon", 
Frankford  Arsenal  Report  M73-30-1,  Oct  1973,  pp  1-32. 
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These  factors  ere  corrected  by  extending  en  analysis  due  Initially 
to  Bland  ^ to  the  present  cartridge  case  problem.  In  Bland's  analysis^ 
the  state  of  stress  and  strain  in  an  elastoplastic,  workhardening 
cylinder  subject  to  internal  and  external  pressure  and  temperature 
was  given.  Because  the  piece-wise  linear  Tresca's  yield  criterion 
was  adopted,  the  resulting  equations  were  linear,  and  solvable  in 
closed  form  for  linear  work  hardening.  Bland  assumed,  however, 
that  the  cylinder  wall  was  so  thick  that,  at  all  levels  of  loading,  part 
of  the  cylinder  was  elastic,  and  part  plastic.  His  formulation  can  not 
readily  handle  the  more  difficult  case  in  which  the  entire  cylinder 
goes  plastic  before  loading  is  complete. 

In  the  prosent  analysis,  it  is  assumed  that  Bland's  analysis  is 
applicable  until  the  entire  case  has  yielded.  Then  Bland's  equations 
are  reformulated  under  the  assumption  that  the  entire  case  ex- 
periences some  plastic  deformation  and,  thus,  that  no  region  remains 
purely  elastic.  Expressions  for  the  states  of  stress,  strain  and  dis- 
placement as  functions  of  internal  pressure  are  presented.  These  ex- 
pressions are  valid  during  the  loading  of  the  case.  Unloading  and  the 
associated  evaluation  of  the  radial  stress  at  extraction  have  not  been 
treated. 


THEORY 

Cartridge  Case  Model 

Consider  the  cartridge  case  to  be  a cylinder  undergoing  axi- 
symmetrlc  loading.  Then  the  radial  equation  of  equilibrium  is 

^ a Q 

dP 


2 

D.  R.  Bland,  "Elastoplastic  Thick- Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temperature 
Gradients",  J.  Mech.  Phys.  of  Solids,  Vol  4,  1956,  pp  209-229. 
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The  email  strain  - small  displacement  equations  are 


• 

t?  - dii 

(2) 

dP 

e 

s VA/  p 

(3) 

where  G(«  0$  VJk.  Bp  B>0*  are*  respectively,  the  radial  and 

tangential  stresses,  the  radial  displacement,  the  radial  and  tangen- 
tial strains. 


The  total  strains  are  assumed  to  be  made  up  of  elastic  and  plastic 
parts,  i.  e.  t 


€ ^ ? 


6?  = 

(4.  a) 

6&  - 

_ e p 

to  ■**  6*6 

(4.  b) 

64  * 6* 

(4.c) 

i 


where  the  superscripts  "e"  and  "p"  denote  elastic  and  plastic  com- 
ponents, and  is  the  total  axial  strain. 

The  elastic  strains  are  related  to  the  stresses  through  three 
Hooke's  Law3  equations  of  the  form 

E.  fcp  = <TP-  + E <*.T  (s.») 

E.  C^-  + * E«*.T  (5.  b) 

E- %=  a#-  wCja'-t-cp  +•  e t (5.0 

3 

A.  E.  H.  Love,  Mathematical  Theory  of  Elasticity,  Fourth  Edition, 
Cambridge  University  Press,  1927,  pg  108. 
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where  E la  the  Young's  modulus  of  th«  material  and  ie  the  Poisson's 

ratio,  both  of  which  are  temperature  insensitive.  T Is  the  temper- 
ature rise  and  Gk  the  coefficient  of  thermal  expansion. 

In  order  to  relate  the  plastic  strain  components  of  Equations  4 to 
the  stresses,  the  following  assumptions  are  made: 

1.  The  material  is  linear  work- hardening  after  yielding. 

2.  The  Tresca  yield  criterion  is  used. 

3.  The  stresses  are  such  that: 

Ge  ^ CTfc  ?/  Op  (6) 


4.  The  difference  between  any  two  stresses  never  drops  below 
the  yield  strength  of  the  material. 

As  a result  of  these  assumptions  Bland2  shows  that  the  associated 
(Tresca)  flow  rule  provides  the  following  expressions  for  the  plastic 
constitutive  equations: 


CVG-e  = K6(.^7&t) 

(7) 

CP  =■  v.2/v3)Lq 

(8) 

& a = 0 

(9) 

% ■=-%  f/P 

(10) 

1 

D.  R.  Bland,  "Elastoplastic  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients",  J.  Mech.  Phys.  of  Solids,  Vol  4,  1956,  pp  209-229. 


where  ^ is  the  yield  stress  in  shear,Ct\fp  i*  the  slope  of  the  plastic 
portion  of  the  stress -strain  curve  in  shear,  is  the  equivalent 
plastic  strain,  and  the  denotes  an  Increment  of  strain. 

An  additional  equation  is  reauired  in  order  to  maintain  equilibrium 
in  the  longitudinal  direction,  Ii  P is  the  net  applied  longitudinal  load, 
including  possible  friction  at  r » b due  to  chamber  interference, 
then 

P =.  at?  I PCTjAr  (n> 

Jo. 


is  the  required  expression  of  longitudinal  equilibrium.  Here  b and  a 
* re,  respectively,  the  outer  and  inner  case  radii. 

It  Is  important  to  note  several  facts  at  this  point: 

1.  The  constitutive  Equations  7 to  10  are  linear  while  still  satis- 
fying an  exact  theory  of  plasticity.  This  is  a considerable  simplifi- 
cation over  the  classical  treatment4  which  used  the  non-linear  von- 
Mlsos  criterion,  and  which  may  be  solved  only  by  a numerical  ( finite  - 
difference)  method. 

2.  The  equilibrium  equations  above,  1 and  2,  cannot  be  solved  ex- 
plicitly to  obtain  the  stresses  as  a function  only  of  the  applied  loads 
without  regard  to  the  constitutive  equations  and  material  properties  as 
well  as  was  done  in  the  previous  case  extraction  study. 4 

3.  Lastly,  because  Equations  2 and  3 are  to  be  solved  simultane- 
ously, the  resulting  total  strains  and  will  be  compatible. 
Physically  this  m-»ans  that  the  cartridge  case  will  now  "fit  together" 
after  deform  ation. 


R*  Hill,  Mathematical  Theory  of  Plastlcit 
pp  109-120. 


Clarendon  Press,  1950, 
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At  this  point,  there  are  14  Equations  (1  to  5,  7 to  11)  in  as  many 
unknowns:  (3  stresses,  3 total  strains,  3 plastic  strains,  3 elastic 
strains,  1 equivalent  strain,  and  1 displacement).  Thus  the  problem 
is  determinate. 

By  combining  Equations  1,  7 and  8 we  know  that 


A value  for  t>b  U given  by  Bland's  ^Equation  (19).  This  when 
combined  with  Equation  12  leads  to  the  following  differential  equation 

for  C(»  s 

p ds  0+  *0 72/tf.  ^‘)  = Ko  V 
|ftoTZ/v3]  |^C/P1v<|5  juc)rai'-TCt)><<U^J  (1„ 

C is,  as  yet,  an  undetermined  constant,  and  T ( r ) is  the  radial  temper- 
ature distribution. 

It  is  assumed  that  the  case  is  heated  non-uniformally  to  a temper- 
ature Ta  (above  ambient)  due  to  the  propellant  at  r = a,  and  is  heated 
to  T^  at  r = b at  the  barrel  case  interface.  The  steady  state  temper- 
ature through  the  wall  is  then  a solution  of  Laplace's  equation,  which 

is 


Ttt)  = Mf  N lap 


(14. a) 


D.  R.  Bland,  "Elastoplastic  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients  ",  J,  Mech.  Fhys.  of  Solids,  Vol  4,  1956,  pp  209-229, 
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where 


M - [_"to.  'at  - "^bWvcJJ  \n.  (&/b) 

N*  [Tb-T«.]  \a(e/b) 


Introduction  o£  Equation  14  into  13  and  integration  leads  to  the 
following  expression  for 


Op- 


(14.  b) 


(14.  c) 


where 


-tf,  - (hs?  2/vS)  ['  + Ct  -V  ‘->J 


(16.  a) 


(16.  b) 


are  material  constants  introduced  for  convenience.  C and  O are 
constants  to  be  determined  by  the  boundary  conditions. 

The  hoop  stress*  £-fe  , is  obtained  from  Equation  7 as 

Ge-  Gr  + ht[i  * 
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The  value  for  C5p  in  Equation  17  is  obtained  from  Bland'  a*  Equation 
19  with  the  aid  of  15.  Substituting  that  result  in  17  the  final  expression 

{or  <5q  u 


Oe=  kct  aK  7^/p2 


* •CCw**)  n]  t ffr 
2 


where  is  given  by  15. 

The  expression  for  the  axial  stress,  , may  be  obtained  from 
Equation  15  and  Bland' s * Equation  13  as 

- E*<T 

where  the  longitudinal  strain,  &£*  is  given  by  ^ 

~[p/ir  +aEoCj  oT<1J  [E  it'-of)] 

where  p and  q are,  respectively,  the  Internal  (propellant)  and  ex* 
ternal  (interference  between  case  and  chamber)  pressures,  or 

p . . <rpcO  , % = - G.-Cb ) 


D.  R.  Bland,  "Elastoplastic  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients",  J.  Mech.  Fhys.  of  Solids,  Vol  4,  1956,  pp  209*229. 
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A*  mentioned  previously,  P is  the  net  axial  force  to  which  the 
cartridge  case  is  subjected  at  extraction.  If  experimental  values  for 
the  thrust  from  the  gas  expulsion  at  the  case  mouth  were  known,  a 
value  of  P could  be  determined.  This  would  in  turn  allow  the  average 
longitudinal  strain  and  stress  i t be  determined.  Such  experimental 
data  is  not  as  yet  available.  * 

A unique  feature  of  this  analysis  is  that  th<*  formula  for  Gr  . 
Equation  IS,  does  not  depend  on  , and,  hence,  does  not  depend 
on  the  axial  thrust,  P.  This  is  a great  simplification  over  the  classic 
treatment4,  in  which  the  value  of  CSV  was  nonlinearly  coupled  to  both 
the  axial  strain  and  axial  stress.  Thus  the  lack  of  data  for  P does  not, 
in  the  present  analysis,  explicitly  influence  the  determination  of  the 
force  of  extraction,  since  the  force  is  a function  only  of  Or  * 4 As 
will  be  shown  below,  however,  the  radial  displacement,  u,  does 
explicitly  depend  upon  Ul  and  P.  During  extraction  u is  used  to 
determine  the  point  of  unloading  and  thus  affects  QV  * 

In  order  to  calculate  u,  the  radial  displacement,  use  is  made  of 
Equations  14  and  15  by  substitution  into  Bland' s Equation  16. 2 
That  result  is; 

- vezr  * 

tO**’}*  [_'v'p  ^ (2] 

The  results  maybe  summarised  as  follows;  The  stresses  are 
given  by  Equations  15,  18,  and  19,  respectively.  The  axial  strain, 

# is  given  by  20,  and  the  radial  displacement,  u,by  Equation  21. 
With  u known  the  total  strains,  , fc,  ($  , follow  from  Equation*  2 
and  1;  the  elastic  strains  C L , follow  from  5.  a and  5.b  and  the 

v j> 


D.  R.  Bland,  "Elastoplastic  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients".  J.  Mech.  Phys.  of  Solids,  Vol  4,  1956,  pp  209-229. 

4R.  Hill,  Mathematical  Theory  of  Plasticity.  Clarendon  Press,  1950, 
pp  109-120. 

5S.  Goldstein,  Private  Communication. 
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f f 

plastic  strains,  from  Equations  4.a-4.c. 

The  above  equations  for  displacements,  strains  and  stresses  con- 
tain two  arbitrary  constants,  C and  D.  As  previously  mentioned  these 
are  determined  by  the  boundary  conditions.  The  boundary  conditions 

are: 

a.  During  expansion  but  prior  to  chamber  cosdaet: 

Of  tcO=-P  <7fCb)so 

b.  During  expansion  and  chamber  contact  but  prior  to  unloading: 


GpCO  - -p 

(23) 

Of  to  - Cfic) 

(24) 

Ul  (.<-)  = U.V) 

(25) 

Of  <.d)-C 

(26) 

where  c and  d are  the  inner  and  outer  chamber  radii,  respectively, 
and  the  primed  (")  quantities  pertain  to  the  revelant  constants  within 
the  chamber. 

For  case  (a),  the  conditions  given  in  Equation  22  and  Equation  15  to 
evaluate  C and  D.  Those  equations  are  in  matrix  form,  [E][F]  - [G], 
where  the  elements  of  the  column  matrix  [F]  are  the  unknowns  C and 
D: 


\ 

C 

-p 

\ 

0 

o 

(27) 


Hence  the  formal  solution  for  F is 

M -•  MM 


(28) 


For  esse  (b),  the  radial  stress  sad  displacement  given  in 
Equations  15  and  22  must  be  used  to  match  the  radial  stress  and 
displacement  in  the  elastic  chamber.  Expressions  for  the  state  of 
stress  and  displacement  in  the  chamber  are  summarised  in  the 
Appendix.  Using  Equations  A-l,  A-4,  15  and  22  yields  in  the  same 
matrix  notation  [£][F]  * [G]t 


1 

O 0 

c 

-^pT^jida) 

o O 

1 -f* 

0 

O 

-» 

1 -t>* 

E 

f 

&4 

where 

^V.i-a*OC bj  E. 

- v'-  ao*/]  /E#  (10) 

E*  * -(.ifvO/e'b 

a ^O-lwXrtv^taTIjJflb/E  “ ■* 

♦ «.Utv)[nt>  rNtblnb-b«)]\ 


I 


I 

l 

M 

■i 


4 

*> 

i 
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The  solution  matrix,  [Fj,  te  again  given  by  Equation  28. 


The  actual  calculating  procedure  it  summarized  as  follows: 

If  chamber  contact  has  not  been  made.  Equation  27  is  used  to  calcu- 
late C and  D.  If  contact  has  been  made.  Equation  29  is  used  to  calcu- 
late C,  D,  K|,  Kg.  The  remaining  variables  are  calculated  as 


Cartridge 

Case 


Variable  From  Equation 


Co 

fei 

cr* 

u. 

fcf 

te 


15 

18 

20 

19 

21 

2 

3 


Chamber 


<j/  A-l 

Oe'  A'2 

Oi 

. > t A-4 


The  values  of  Oq  , , and  Or  calculated  above  are  used  to 

check  the  inequality  given  by  Equation  6 and  also  that  no  difference 
between  any  two  of  the  stresses  drops  below  the  yield  stress  . 
A discussion  of  these  requirements  ia>  given  in  Bland  . If  either  of 
these  requirements  is  not  fulfilled  the  analysis  is  not  valid. 


As  mentioned  in  the  Introduction  the  present  analysis  is  valid 
only  during  loading  of  the  case,  unloading  has  not  as  yet  been  treated 
for  this  formulation.  Unloading  for  crude  theories  of  plasticity  has 


D.  R.  Bland,  "Elastoplastic  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients",  J.  Mech,  Phys,  of  Solids,  Vol  4,  1956,  pp 209-229. 


been  treated  approximately  in  other  studies 
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Transient  Heating 

The  steady  state  heat  transfer  problem  is  presented  above.  A 
more  accurate  solution  can  be  obtained,  however,  by  considering  the 
propellant  heating  to  be  transient. 

The  transient  temperature  distribution,  T(r,t),  is  assumed  to  be 
known  by  separate  analysis.  References  8 and  9 give  the  details  for 
such  calculations.  Equations  13  may  be  integrated  in  general  formtoyield. 

CpCfyO*  19.0^  )kr  - *iC/arJ  0 + 

K.w.utvij  (J1> 


for  the  traneient  radial  stress.  Expressions  for  CTftiUjblit  0 
follow  closely  the  analysis  given  in  the  proceeding  section. 

RECOMMENDATIONS 


A theoretical  analysis  of  the  states  of  stress,  strain  and  displacement 
n.  * cartridge  case  and  chamber  has  been  performed.  The  cartridge 
case  has  been  modeled  as  a right  circular  cylinder  of  thermo  elastic- 
plastic,  linear  strain  hardening  metal.  The  chamber  is  elastic. 


V.  Gordon,  "Analytic  Study  of  Extraction  In  the  M16  Weapon", 
Frankford  Arsenal  Report  M73-30-1,  Oct  1973,  pp  1-32. 

^L.  M.  Gold,  "Cartridge  Case  Chamber  Interaction  During  Firing", 
Frankford  Arsenal  Report  M73-35-1,  December  1973. 

®A.  Carslaw,  J,  Jaegar,  Conduction  of  Heat  in  Solids,  Second  Edition 

Oxford  University  Press,  195$, 
a 

B.  Boley  It  J.  Weiner,  Theory  of  Thermal  Stresses,  New  York, 

John  Wiley,  i960. 
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It  was  found  that  by  extending  an  analysis  of  Bland2  the  short- 
comings In  a previous  simplified  model*  could  be  corrected.  These 
corrections  included! 

1.  Thermal  strains  due  to  propellant  heating. 

2.  Compatibility  of  the  radial  and  hoop  strains. 

3.  Improved  representation  of  plastic  flow  (consltutlve  equations). 

It  was  also  noted  that  the  present  analysis  Is  still  far  simpler 
than  the  classical4  analysis,  which  requires  extensive  numerical  In- 
tegration of  nonlinear  equations.  The  present  analysis  reduces  to 
algebraic  expressions  for  the  stresses,  strains  and  displacements, 
requiring  only  the  Inversion  of  either  a2x2ora4x4  matrix. 

While  the  formulation  has  been  specifically  used  for  a steady  state 
temperature  effect  (via.  Equation  14),  the  extension  to  transient 
heating  is  demonstrated. 

It  is  recommended  that: 

1.  This  analysis  be  programmed  and  the  results  compared  with 
experiment  when  available. 

2.  The  analysis  be  extended  to  include  unloading  and  the  actual 
calculation  of  extraction  forces. 


V.  Gordon,  "Analytic  Study  of  Extraction  in  the  M16  Weapon", 
Frankford  Arsenal  Report  M73-30-1,  Oct  1973,  pp  1-32. 

ZD.  R.  Bland,  "Elastoplastlc  Thick-Walled  Tubes  of  Workhardening 
Material  Subject  to  Internal  and  External  Pressures  and  Temper- 
ature Gradients",  J.  Mach.  Phys,  of  Solids,  Vol4,  1956,  pp  209-229. 

*R.  Hill,  Mathematical  Theory  of  Plasticity.  Clarendon  Press,  1950, 
pp  109-1*5: 
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APPENDIX 


State  of  Stress  end  Displacement  In  the  Chamber 


Consider  the  chamber  to  be  a long  right  circular  annulus, 
c £ r 4k  d,  in  a state  of  plane  strain.  Than  the  most  general 
solutions  to  the  elastic  isothermal  governing  equations*  are 

Of  = k/-  w't  / p1 

(A-l) 

Oils  * Kf  «■ 

(A-2) 

C V = a^'K,' 

(A- 3) 

s 'V'f’O'^’-at.voO/e 

+ r (a-4) 

where  Kj  and  K2  are  consorts  to  be  determined. 

If  the  chamber  is  too  hot  for  the  tixermcl  strains  in  the  chamber 
to  be  neglected  than  A-l  through  A-4  can  be  easily  aiollfled  for 
thermal  strain  as  shown  by  Medal. 7 


7 

A.  Nadal,  Theory  of  Flow  and  Fracture  in  Solids,  Vol  2,  McGraw- 
Hill,  1963,  pg  389. 
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